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Abstract. D0 provides a weahh of measurements conceived for probing perturbative and 
non-perturbative aspects of QCD, giving an accurate experimental account for Standard Model 
production processes including jets, leptons and photons and improving the sensitivity and the 
understanding in searches for new physics. Among the most important subjects are inclusive 
jet production, vector boson plus jet production, direct photon production and measurements 
in Minimum Bias events and Double Parton Scattering which are discussed here and compared 
to theory predictions. 



1. Introduction 

Proton anti-proton collisions at a centre of mass energy of ^/s = 1.96 TeV are taking place at the 
the Tevatron accelerator, Fermilab. A peak luminosity of 4.3 • 10^^ cm~^s~^ has been reached 
and over 11.6 fb^^ integrated luminosity are already delivered, of which 10.5 fb^^ have been 
recorded by the D0 experiment. The presented analyses are based on integrated luminosities 
ranging between 0.7 and 4.2 fb~^. 

The D0 detector [1] has broad particle identification capabilities. Jets, electrons and photons 
are detected by the calorimeter with fine granularity (in pseudorapidity and azimuthal angle of 
Ar] X A(j) ~ 0.1 X 0.1) and good energy resolution. Charged particles are reconstructed by means 
of the central tracking system and muons are detected in the muon spectrometer as outermost 
detector component. The data taking efficiency is > 90% (in Run 2b it amounts to about 92%, 
on average). 

Jet measurements are well understood by Quantum Chromodynamics (QCD) calculations. 
They can be exploited to set constraints on the non-perturbative Parton Distribution Functions 
(PDF's). In the kinematic plane of proton momentum fraction x and virtuality scale Q^, regions 
accessible to fixed target experiments, HERA, the Tevatron and LHC are complementary. Only 
Tevatron inclusive jet data provide significant constraints at high x and Q^. Jet measurements 
are also important within searches for new phenomena. Heavy particles can decay into quarks 
which in turn fragment into jets. Extra dimensions and quark compositeness are examples of 
such new kinds of phenomena. 

The production of dijets at the Tevatron is dominated by the gluon gluon fusion process 
at low transverse jet momenta. Toward larger transverse jet momenta the fraction of quark 
anti-quark annihilation dominates, while the fraction of quark gluon fusion tends to decrease 
marginally. 

Among all past and present hadron colliders the Tevatron inclusive jet production data 
dominates together with published LHC results the high jet transverse momentum reach. At 



the same time, the Tevatron has by far the best sensitivity at large (anti-)proton momentum 
fractions x. 

While the data are corrected for detector effects, the predictions take into account non- 
perturbative contributions due to the underlying event and hadronisation. Finally, the 
comparison to data is then established at the hadronic final state as defined in [1]. 
The convention to set the speed of light c = 1 is adopted throughout this article. 

The following section covers inclusive jet production measurements. Section [3] shows two 
examples of vector boson plus jet production measurements. Section U] encompasses the 
most relevant direct photon production measurements of D0. Important non-perturbative 
measurements in the regime of minimum bias events and double parton scattering are given 
in section [SJ 

2. Inclusive jet production 

The measurement of the cross section of inclusive jet production in hadron collisions provides 
stringent tests of QCD. Large jet transverse momenta with respect to the beam axis (px) ensure 
small contributions from long-distance processes and the jet production can be calculated in 
perturbative QCD (pQCD). The inclusive jet cross section in pp collisions provides one of the 
most direct probes of physics at small distances. In particular, it is directly sensitive to the 
strong coupling constant (a^) and the PDF's of the proton. Furthermore it can be exploited to 
set constraints on the internal structure of quarks [5] . 

2.1. Inclusive jet production 

The double differential inclusive jet production cross section is measured [2] as a function of the 
leading jet transverse momentum and the rapidity based on an integrated luminosity of 0.7 fb^^. 
The D0 Run II iterative seed-based cone jet algorithm including mid-points [3] with cone radius 
R = \J (Ay)'^ + (A(/))2 = 0.7 in rapidity y and azimuthal angle (p space is used to cluster energies 
deposited in calorimeter towers as well as for the pQCD calculations on the partonic final state. 
The jet transverse momentum is corrected for the energy response of the calorimeter, energy 
showering in and out of the jet cone and for additional energy from event pile-up and multiple 
parton interactions. Events with jets above a transverse momentum threshold of 50 GeV in the 
rapidity range of \y\ < 2.4 are selected. The cross section is measured in six rapidity bins as a 
function of jet pT, reaching transverse momenta above 600 GeV. The cross section is corrected 
for muons and neutrinos, not reconstructed within jets. Corrections for jet migration between 
bins in pT and y due to finite resolution in energy and position are determined in an unfolding 
procedure, based on the experimental pT and y resolutions. The bin sizes are chosen to minimise 
migration corrections due to the experimental resolution. The dominant systematic uncertainty 
comes from the Jet Energy Scale (JES). The steeply falling jet transverse momentum spectrum 
causes shifts in energy which translate into a large uncertainty on the cross section. The typical 
JES uncertainty is 1 - 2% which results in a total uncertainty on the cross section of 15 - 
30%. Perturbative QCD predictions to next-to- leading order (NLO) in as, computed using the 
fastNLO program [6j, based on NLOJET++ [7] are compared to data making use of PDF's 
from CTEQ6.5M [8j |9]|10). The theoretical uncertainty due to variation of the renormalisation 
(jiR) and factorisation (hf) scales by a factor of two up and down are of order 10%. The scale 
itself is chosen to be the individual jet pT- In all y regions the predictions agree well with the 
data. The data tends to be lower than the central prediction making use of CTEQ PDF's, in 
particular at very large transverse jet momentum. 

2.2. Inclusive dijet production invariant mass cross section 

The inclusive dijet production cross section is measured [11] as a function of the dijet mass 
based on an integrated luminosity of 0.7 fb~^. This measurement is sensitive to quark 



compositeness, extra spatial dimensions and undiscovered heavy particles decaying into two 
quarks. The invariant mass distribution is particularly sensitive to the PDF of gluons at high 
proton momentum fraction x where the gluon distribution is weakly constrained. The two 
leading jets (cone jet algorithm [3] with radius R = 0.7) are required to exceed a transverse 
momentum threshold of 40 GeV. The cross section is measured in different regions of maximal 
rapidity, defined as |y|max = max(|yi|, \y2\)- The six equidistant adjacent regions are given by: 
< |y|max < 0.4, 0.4 < lyUax < 0.8, 0.8 < lyUax < 1-2, 1.2 < lyUax < 1-6, 1.6 < |y|max < 2.0 
and 2.0 < \y\ 

max < 2.4. Calorimeter shower shapes are used to remove remaining background 
due to electrons, photons and detector noise giving rise to jets. The invariant mass bins are 
chosen to have a size of about twice the mass resolution and to respect an efficiency and purity 
of about 50%. The jet pr resolution is measured in events with exactly two jets and has been 
found to be approximately 13% (7%) at pT — 50 (400) GeV in the Central Calorimeter (CC) 
and the Endcap Calorimeter (EC) and 16% (11%) at pT — 50 (400) GeV in the Inter Cryostat 
Region (ICR) between the two calorimeter compartments CC and EC. The total experimental 
corrections to the data are less than ±2%. The systematic uncertainties on the cross section are 
dominated by the uncertainties in the jet energy calibration which range from 6 - 22% in the EC 
to 15 - 45% in the CC. The second largest systematic uncertainty is given by the px resolution 
uncertainty, ranging between 2 and 10% in all regions. The luminosity determination has an 
uncertainty of 6.1% which is completely correlated across all bins. The data are compared to the 
NLO prediction computed using [6], based on NLOJET++ [7|[l3] for MSTW2008NLO PDF's 
[12J with as{mz) = 0.120. In addition, a comparison to theoretical predictions using CTEQ6.6 
PDF's [14J has been made. The difference in cross section due to the choice of PDF's is 40 - 
60% at the highest mass. MSTW2008NLO PDF's are favoured. However, it has to be noted 
that their determination included the D0 inclusive jet production cross section measurement [2] 
which is based on the same data set as this dijet measurement. A further difference is that the 
MSTW2008NLO PDF's exclude Tevatron data taken before 2000 in contrast to the CTEQ6.6 
PDF's. 



2.3. Inclusive 3- jet invariant mass cross section 

The inclusive 3-jet production cross section is measured [IS] as a function of the invariant 3-jet 
mass based on an integrated luminosity of 0.7 fb~^. Recent measurements of inclusive jet and 
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Figure 1. Data over theory ratios of the differential cross section da^jet / dM^j^t in the different 
rapidity regions and for different third leading jet transverse momentum requirements. The 
pQCD predictions are given for two different PDF's. 



dijet production in pp collisions at a centre of mass energy of t/s = 1.96 TeV [16] [T7] [18] [2] [19] [20] 
have been used to determine ag [21] and the proton PDF's and to set limits on a number 
of physics models beyond the standard model |22j |23) (23]. This shows the successful pQCD 
description of observables which are directly sensitive to matrix elements of 0{al). Testing 
pQCD at higher orders of as requires measuring cross sections of higher jet multiplicities. The 
three jet cross section is directly sensitive to matrix elements of ) while the PDF sensitivity 
is similar. Due to the fact that the pQCD calculations are available to next-to-leading order 
(NLO) in Us [25] |2S] [Z] [E] ) the three jet cross section can be used for precision phenomenology, 
e.g. the simultaneous determination of ag and PDF's from experimental data. Together with 
the inclusive jet and dijet cross sections the results for as and the PDF's get then partially 
decorrelated. 

The leading transverse momentum jet has to fulfil px > 150 GeV and the third leading 
transverse momentum jet px > 40 GeV. The leading jets are restricted to \y\ < 0.8, \y\ < 1.6, 
\y\ < 2.4 in three different measurements. Two additional measurements are made for a third 
leading jet pT > 70 GeV and pr > 100 GeV in the rapidity range of \y\ < 2.4. The cone jet 
algorithm [3] with radius R = 0.7 is used and all pairs of the three leading jets are required 
to be spatially separated by AR > 1.4(= 2 ■ R). Dominating systematic uncertainty on the 
cross section is given by the Jet Energy Scale (JES) with 10 - 15% for \y\ < 0.8 and 10 - 
30% for |y| < 2.4. The px resolution for jets is about 15% at 40 GeV, decreasing to below 
10% at 400 GeV. The corresponding uncertainty on the cross section amounts to 1 - 5%. The 
uncertainty on the luminosity contributes with 6.1%. In Fig. [1] the data over theory prediction 
ratios of the differential cross section da^^ct /dM^jet are shown in the different rapidity regions and 
for different third jet transverse momentum thresholds. The MSTW2008NLO PDF uncertainties 
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Figure 2. Trijet inclusive over dijet inclusive cross section ratio as a function of the leading 
transverse momentum jet in comparison to various predictions of SHERPA and PYTHIA (top) 
and next-to-leading order pQCD calculations (bottom). 



corresponding to the 68% C.L. are shown by the hght band. Normahsed to MSTW2008NLO, 
the CTIO PDF's[23] with corresponding value of as{Mz) = 0.118 predict a different shape with 
increasing discrepancy at higher invariant masses. 

2.4- Ratio of inclusive trijet/dijet production 

The ratio of inclusive trijet over dijet production is measured [28] as a function of the leading 
jet transverse momentum p™*^^ based on an integrated luminosity of 0.7 fb^^. Jet production is 
sensitive to the dynamics of the fundamental interaction and the partonic structure of the initial 
state hadrons. Measurements dedicated to investigate the interaction dynamics are preferentially 
based on observables which are insensitive to the PDF's. Such observables can be constructed 
as ratios of cross sections for which the PDF sensitivity cancels out. The ratio i?3/2 of inclusive 
trijet over inclusive dijet production is such an observable. In the R3/2 measurement reported 
here is being probed up to jet transverse momenta of 500 GeV. Jets are defined by the cone 
jet algorithm |3j with radius R = 0.7 and all pairs of n-leading jets (n = 2, 3) are required to be 
separated by AR > 2 ■ R. The cross section ratio is determined in three bins of jet transverse 
momentum thresholds > 50, 70, 90 GeV. The leading transverse momentum jet has to 
exceed the n-th leading jet by p™'^^ > p™™ + 30 GeV to ensure sufficient phase space for the 
second and third leading jet such that corrections due to the experimental pT resolution remain 
small. In Fig. [2] (top) the cross section ratio is compared to the event generators SHERPA 
1.1.3 [29J and PYTHIA 6.419 [27| for three different tunes using a pT-ordered parton shower. 
The MSTW2008LO PDF's have been deployed for the predictions of both event generators. 
While SHERPA with leading order matrix elements for 2-, 3- and 4-jet production matched 
to parton shower describes the data well, PYTHIA with tunes using the py-ordered parton 
shower is not able to describe the data. PYTHIA making use of the tune BW with Q^-ordered 
parton shower is able to describe the data but a former D0 measurement of dijet azimuthal 
decorrelations [30] can not be described by the BW tune. As can be seen in Fig. [2] (bottom) 
the prediction of pQCD in next-to-leading order is describing the data very well. PYTHIA 
predicts non-perturbative corrections due to hadronisation and underlying event of 2.5% it 
0.5%. Therefore the discrepancies of various tunes can not be explained by poor parameter 
choices for the hadronisation and/or underlying event models. It can only be explained by 
parameters which affect the perturbative physics as implemented in the parton shower or by a 
fundamental limitation of the model itself. 

3. Vector boson plus jet production 

Measurements of vector boson plus jet production provide fundamental tests of pQCD in 
particular at high momentum scales. Furthermore these kind of processes can constitute the 
dominant background in measurements of single top quark and ti production as well as in 
searches for the Standard Model (SM) Higgs boson and for physics beyond the SM. In particular 
Z boson plus b jet production is an important background for associated Standard Model Higgs 
boson production {ZH — )• Zbb) and also super symmetric partners of b quarks. Furthermore the 
b quark density of the proton can be probed. 

3.1. Inclusive Z boson cross section ratio cj{pp Z + b + X)/a{pp — )• Z + j + X) 
The ratio of the inclusive production of a Z boson accompanied by a 6 jet over the inclusive 
production of a Z boson accompanied by a light flavoured jet a{pp Z + b + X)/a{pp — )• 
Z + j + X) is measured [31] as a function of the leading jet transverse momentum based on 
an integrated luminosity of 4.2 fb~^. The measurement of the cross section ratio benefits from 
cancellations of many systematic uncertainties. 

Selected events are required to pass at least one of the single electron or single muon triggers. 
A Z boson candidate with an invariant dilepton mass of 70 < mu < 110 GeV has to be 



reconstructed. Electrons (muons) have to have a transverse momentum of px > 15 GeV {pT > 
10 GeV). Electron candidates are built from isolated energy deposits in the electromagnetic 
calorimeter with a shower shape consistent to that expected from electrons. Combined tracking 
and calorimeter isolation requirements are applied to the muon candidates. At least one jet 
with pt > 20 GeV and further jets with px > 15 GeV, making use of the cone jet algorithm [3] 
with cone radius R = 0.5 are required. Jets with b content have a different energy response and 
receive an additional average energy correction of about 6%. The 6-tagging algorithm is based on 
a neural network exploiting the longer lifetime of 6-flavoured hadrons which results in tracks with 
sizeable impact parameters and secondary vertices. The necessary tracks associated to jets have 
to exceed a transverse momentum of 0.5 GeV and 1.0 GeV for the leading track. At least one of 
the jets has to be 6-tagged. The tagging efficiency for b jets and light jets are parameterised as 
functions of jet pT and pseudorapidity rj. They amount to about 58% and 2% respectively. To 
further separate b jets from c and light jets a discriminant based on the secondary vertex mass 
Msv and the jet lifetime impact parameter (JLIP) is introduced. Events with missing transverse 
energy > 60 GeV are rejected to suppress background from production. The dominant 
background to Z + jet production arises from multijet events with jets misreconstructed as 
leptons, especially in the ee channel. Smaller background contributions arise from and diboson 
{ZW, ZZ, WW) production. A binned maximum likelihood fit to the discriminant distribution 
in data using a linear combination of light, c and b flavoured jet templates from simulation is 
performed to determine the flavour fractions. Several experimental uncertainties cancel out in 
the measurement of the cross section ratio a{Z +b+ X) / a{Z + j + X) , including the uncertainties 
on the luminosity, trigger, lepton identification and some jet identification efficiencies. The two 
largest remaining sources of systematic uncertainty are due to the discriminant efficiency and the 
shape of the discriminant templates. Other important sources of uncertainty are the 6-tagging 
efficiency with 2.4%, the 6 jet energy scale with 2% and reconstruction efficiency with 3.2%. The 
measurement of the ratio yields ^j^'^jg^^^ = 0.0193 ± 0.0022(stat) ± 0.0015(syst) consistent with 
the ratios obtained separately in the two lepton channels and the NLO prediction of MCFM |32] 
of 0.0192 lb 0.0022 for renormalisation and factorisation scales fij^ = Hp = making use of the 
MSTW2008 [12J PDF's. 

3.2. Inclusive W{-^ eVe) plus jets production rates 

The inclusive W boson plus jets production rates are measured [3B] in dependence of the inclusive 
number of jets based on an integrated luminosity of 4.2 fb~^. The results are corrected to the 
hadronic final state in deconvolving the effects of finite detector resolution, detector response, 
acceptance and efficiencies. Jets are identified with the D0 midpoint cone algorithm [3] with 
cone radius R = 0.5. Jet energy corrections are determined by means of the transverse 
momentum imbalance in 7 + jet events, where the electromagnetic calorimeter response is 
calibrated using Z/j* — )• e~^e~ events. The missing transverse energy is corrected for the 
presence of any muon in the event. Central electrons with pT > 15 GeV, missing transverse 
energy ^ j- > 20 GeV, transverse mass of the W boson candidate > 40 GeV and jet 

transverse momenta of pT > 20 GeV are required. Furthermore, the distance between electrons 
and the nearest jet has to satisfy AR > 0.5 in azimuthal angle and pseudorapidity r] space. 
After this selection there are backgrounds from Z + jets, tt, diboson, single top quark and 
multijet production. W/Z+iets and tt processes are simulated with ALPGEN v2.11 [38j, 
interfaced to PYTHIA v6.403 [37] for the simulation of initial and final state radiation and for 
the hadronisation. The PYTHIA event generator is used to simulate diboson production. Single 
top quark production is simulated via the COMPHEP [33] generator interfaced to PYTHIA. 
The cross sections for W/Z + jet production used for acceptance and smearing corrections are 
simulated by ALPGEN, corrected with a constant factor to match the inclusive W/Z + jet cross 
section calculated at NLO [M][35]- Additional corrections are applied to W/Z + heavy flavour 
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Figure 3. Ratio of pQCD predictions to measured differential cross sections for the pT of the 
n-th leading jet in I^ + n jet production (left: n = 1 (top), n = 2 (bottom). Right: n = 3 (top), 
n = 4 (bottom)). The shaded areas indicate the theoretical uncertainties due to variations of 
the factorisation and renormalisation scales. 



jet production to match the prediction of NLO calculations. Multijet background is determined 
by means of a data-driven method [39j in which data in a control region orthogonal to the 
data selection is used to determine shape and overall normalisation of the multijet contribution. 
After background subtraction the data are unfolded, including an unfolding bias correction of 
the order 0.5 - 2%. All differential cross section measurements are normalised to the measured 
inclusive W boson cross section. The total inclusive W boson cross section is measured to be 
aw = 1132 lb l(stat)]l^g4(syst) + 69(lumi) pb with W — )• — ^ efeZ^r contributing as signal. 
A recent publication ^40j of + 4 jet production at NLO for pp collisions at ^/s = 7 TeV is 
available. These predictions are not available for the Tevatron. Comparisons with theory are 
therefore limited to LO for W + 4 jet production. For the predictions up to three partons in 
the final state the partonic final state generators BLACKHAT [41], interfaced to SHERPA [29] 
and ROCKET [42] interfaced to MCFM [32j are used. BLACKHAT + SHERPA deploys the 
renormalisation and factorisation scale fi = fXR = fip = \Ht where Ht is the scalar sum of 
parton and lepton transverse energies. The choice made by the ROCKET + MCFM authors is 

^ = \J M^r + \ {^P^^^Y in the n = 2, 3 and 4 jet bins. Li the 1 jet bin the factor 1/4 is omitted. 
The PDF's used for the theory calculations are MSTW2008 p!2|. Hadronisation corrections are 
obtained with SHERPA and the CTEQ6.6 PDF's. Fig. [3] shows ratios of the pQCD predictions 
to the differential cross sections for W boson plus jet production in the jet multiplicity bins 
n = 1 — 4 as a function of the n-th jet transverse momentum respectively. The BLACKHAT 
-|- SHERPA predictions are in very good agreement with the data everywhere except at small 
jet pt- The ROCKET -|- MCFM predictions are significantly below the data. In particular the 
difference in the W + 2 jet bin indicates that the scale uncertainties are larger than determined. 
The data offers here guidance for good choices of the renormalisation and factorisation scales. 

4. Direct photon production 

Direct photons come unaltered from the hard subprocess and provide therefore a direct probe 
of the hard scattering dynamics. The energy calibration of photons is better than the energy 
calibration of jets. Direct photons are also important background to many physics processes. At 
the same time the cross section of multijet production dominates. Therefore the understanding 
of the QCD production mechanism is indispensable in search for new physics. 



Photon candidates are given by electromagnetic showers in the calorimeter. By means of 
shower profile and geometric isolation criteria they are distinguished from background, such as 
neutral mesons and jets with high electromagnetic fraction. No associated track has to be found 
and jets have to be separated by a geometrical distance of AR > 0.7 in (rj, (/))-space. Multiple 
preshower hits help to distinguish prompt photons from vr*^ — t- 27 decays. 

At low photon transverse momentum (pj, < 120 GeV) photon + jet production through 
Compton scattering qg — )• dominates. The requirement of isolated photons reduces the 
contribution of photon fragmentation of quarks. Photon + jet production can be exploited to 
probe PDF's at low x (down to x ~ 0.007 at the Tevatron), where quarks are constrained by 
HERA data and therefore one is sensitive to the gluon density. Finally pQCD can be probed at 
NLO, including soft gluon resummation and models of gluon radiation. 



4-1- Inclusive photon plus jet production 

The inclusive photon plus jet production is measured [33] as a function of the photon transverse 
momentum based on an integrated luminosity of 1.0 fb~^. An isolated central photon with 
Pt > 30 GeV is required to pass a neural network background discrimination based on tracker 
and calorimeter information. Jets above a transverse momentum of px > 15 GeV in the two 
rapidity intervals |y| < 0.8 (central) and 1.5 < \y\ < 2.5 (forward) are selected. Four regions 
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Figure 4. Ratio of measured triple differential cross section for same/opposite photon and jet 
rapidities and central or forward leading jets as a function of the photon transverse momentum. 



are distinguished, depending on if the photon and jet have the same or opposite rapidity sign 
{y^-y^'^^ < or > 0) and depending on the central or forward location of the leading jet. The triple 
differential cross section as a function of the photon transverse momentum and rapidity and the 
jet rapidity a / dpj.dy'^ dy^^^ is measured. The NLO prediction is given by JETPHOX [4^ [45j 
making use of CTEQ6.5 [8] PDF's and BFG photon fragmentation functions [l6]. In Fig. [3] 
the triple differential cross section is shown for the four different phase space regions, depending 
on the same/opposite photon jet hemisphere and the central/forward leading jet location as a 
function of the photon transverse momentum. The theory is not able to describe the data in the 
whole measured range. The structure is similar to previous observations of UA2 [U], CDF [IS] 
and D0 j49j }50] . The cross sections as a function of the photon pT can also be arranged as 
ratios over two different regions (out of the four possible regions) giving rise to six different 
ratios which benefit from reduced systematics. The shapes of the measured cross section ratios 
in data are qualitatively reproduced by the theory in general, but quantitative disagreement is 
observed for some kinematic regions. This disagreement holds in particular for central jets over 
same sign rapidity forward jets. 

4-2. Inclusive photon plus heavy flavour jet production 

The inclusive photon plus heavy flavour jet production is measured [51j as a function of the 
photon transverse momentum based on an integrated luminosity of 1.0 fb~^. The QCD Compton 
photo production process dominates for photon transverse momenta of px ~ 90 (120) GeV for 
h (c) quark flavour. Since the outgoing quark flavour is for this production mechanism equal to 
the incoming quark flavour, constraints can be set on the heavy flavour content of the PDF's. 
At larger photon pT the quark anti-quark annihilation process with a photon and a gluon 
dominates, where the gluon splitting into a charge conjugated heavy flavour quark pair enters 
into the considered data sample. An isolated central photon with transverse momentum of 
Pt > 30 GeV is required together with central jets (cone radius R = 0.5) above a transverse 
momentum of 15 GeV. A neural network is applied to the jets to enhance the heavy flavour 
content of the selected data sample. Templates for different flavour contents are built by means 
of a function of jet track probabilities. The distributions of these functions are fitted to the 
shape of the data to extract the flavour composition. The measurement is done separately in 
the two kinematic regions given by same and opposite photon and jet rapidities. The cross 
sections for the b and c quark flavours in the two kinematic regions are measured as a function 
of the photon transverse momentum. Comparison to theory [52j reveals disagreement for 7 + c 
jet production for photon transverse momenta above 70 GeV. 

4.3. Direct photon pair production 

The direct photon pair production cross section is measured [53J as function of the invariant 
diphoton mass M^-y, the diphoton transverse momentum pj^, the azimuthal angle between the 
two photons A(j)^^ and the cosine of the diphoton scattering angle in the diphoton rest frame 
cos^* = tanh[(r7i — r/2)/2], based on an integrated luminosity of 4.2 fb~^. Direct photon pair 
production with large invariant mass m^^ constitutes a large and irreducible background to 
searches for the Higgs boson decaying into a pair of photons at hadron colliders. It is also a 
significant background in searches for new phenomena, such as new heavy resonances [54j , extra 
spatial dimensions [55J , or cascade decays of heavy new particles [SGj . Double photon production 
is also interesting to check the validity of pQCD and soft gluon resummation implemented in 
theoretical calculations. Background contributions from single and double photon fragmentation 
are suppressed by photon isolation and invariant m-y-y requirements, which entails smaller 
theoretical uncertainties. Events with two central isolated photons of transverse momentum 
Pj^ > 21 GeV and pj^ > 20 GeV are selected with unbalanced thresholds to ensure non- vanishing 
phase space at NLO. The invariant mass is required to exceed the transverse momentum of the 



diphoton system m^^ > p7p . The two photons have to be separated by the distance Ai? > 0.4 
in (r/, (^)-space. The cross sections are measured differentially as a function of m^^, the diphoton 
transverse momentum ^ the azimuthal angle between the photons Ac/)^^ (see Fig. [5] for a 
data comparison to theory and various models) and the cosine of the scattering angle cos0*, 
where holds p^^ > pj?. The shapes of the pj^ and Acp^^ distributions are mostly affected by 
initial state gluon radiation and fragmentation effects. Furthermore, the m^^ spectrum is in 
particular sensitive to potential contributions from new phenomena. The cos 9* distribution 
probes PDF effects and the angular momentum of the final state. The measured cross sections 
are compared to theoretical predictions from SHERPA [29j, RESBOS [5Z|I5S], DIPHOX [SS] 
and PYTHIA [BD]. Both, RESBOS and DIPHOX provide pQCD NLO predictions, however 
the gg — )• 77 contribution is implemented only at leading order (LO) in DIPHOX. The explicit 
parton-to-photon fragmentation functions are included in DIPHOX at NLO, while in RESBOS a 
function approximating rates from NLO fragmentation diagrams is implemented. Furthermore 
only in RESBOS the effect of soft and collinear initial state gluon emission is resummed to 
all orders. The normalisation uncertainty on the measured cross section amounts to 7.4% 
due to uncertainties on the luminosity of 6.1% and on the selection criteria of 4.3%. PDF 
uncertainties are determined by means of DIPHOX, making use of the CTEQ6.6M [14] PDF's. 
The renormalisation, factorisation and fragmentation scales are set to = fip = fj,f = m^^. 
Studies performed using DIPHOX indicate that the contribution to the overall cross section 
from one- and two-photon fragmentation processes does not exceed 16% and significantly drops 
at large and small A(j)^^ to 1 - 3%. None of the theoretical predictions considered is 

able to describe the data well in all kinematic regions of the four variables. RESBOS shows the 
best agreement with data, although systematic discrepancies are observed at low m^^, high pjJ 
and low A(p^^. The large discrepancy between RESBOS and DIPHOX in some regions of the 
phase space is due to the absence of all-order soft gluon resummation and the LO approximation 
of the gg — )• 77 contribution in DIPHOX. Further insight can be obtained in considering double 
differential cross sections. Therefore the differential cross sections as functions of pj^, A(j)^^ 
and I cos 0*1 are measured in three invariant m^^ mass bins 30 - 50 GeV, 50 - 80 GeV and 80 
- 350 GeV. The largest discrepancies between data and RESBOS are found for each kinematic 
variable in the lowest m^^ < 50 GeV region. This is the region, where the contribution of 



99 ~^ 77 is expected to be largest. Finally, the addition of NNLO corrections to RESBOS 
remains to be compared, once available. 
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Figure 5. Differential diphoton production cross section as a function of the invariant diphoton 
mass m^^ (a), the transverse diphoton momentum pj! (b) and the azimuthal angle between the 
photons A(/)^^ (c) in comparison to theory and various models. 



5. Minimum bias and double parton scattering 

The non-perturbative nature of low energy QCD processes makes it necessary to develop heuristic 
descriptions. Models that describe these soft processes are implemented in Monte Carlo event 
generators. They contain many parameters that need to be tuned using experimental input. 
Most prominent input observables are charged particle multiplicities and track mean transverse 
momenta. 

Double Parton Scattering (DPS) may constitute a background for many rare processes 
of interest, especially those with multijet final states. The measurement of DPS provides 
complementary information about the proton structure in terms of the spatial distributions 
of partons. Furthermore it allows one to probe for parton-parton correlations and therefore to 
check for a possible impact on the PDF's. Multiple Parton Interaction (MPI) models can be 
constrained and improved taking DPS measurements into account. 

5.1. (p and r] correlations in minimum bias events 

The azimuthal angle and pseudorapidity correlations in minimum bias events are measured [61] 
based on a data set taken in the years 2002 to 2006. Complementary to existing observables 
the azimuthal angle difference A(f) between the highest transverse momentum track and any 
other track in a given event is considered here. This observable has the advantage of being as 
independent as possible of the track fake rate and the tracking efficiency. The general shape of 
the A(j) distribution consists of peaks at and vr on top of a pedestal of certain height. 

Events with more than one reconstructed proton anti-proton interaction per bunch crossing 
are selected and the single collision that fired the trigger is being identified. The remaining 
collisions of the same bunch crossing are regarded as minimum bias events. The triggers used 
are dimuon triggers. The primary vertices of these collisions are required to have at least five 
tracks with transverse momentum above 0.5 GeV inside the geometrical tracker acceptance 
(|?7| < 2). Long-lived particles, such as pions and Ks are discriminated by means of a vertex 
requirement based on three dimensional track impact parameter significances. Uncorrelated fake 
tracks and wrongly associated real tracks are distributed uniformly in Acp. This fiat background 
is eliminated in fitting the distribution by a polynomial and subtracting the height at minimum 
from the whole distribution. The absolute level of the shape of the Acp distribution is affected 
by the tracking efficiency but the dependence is minimised by normalising the distribution to 
unit area and only considering the shape. As it turned out it has not been necessary to unfold 
the experimental result. A further observable has been defined by assigning the tracks into 
two pseudorapidity regions based on the sign of the leading tracks r]. The remaining tracks 
belong then to the same or opposite region depending on the sign of their r] respectively. The 
number of same sign tracks subtracted by the number of opposite sign tracks in a given Acj) bin, 
normalised to the sum of all such expressions over all Acj) bins is then the final observable which 
is being considered as a function of A(f). The observables are compared to three PYTHIA tune 
and model implementations, which are Rick Field's Tune A [62], the Perugia tune (PO) |63) 
and the Generalised Area Law model of colour reconnections (GAL) [64j . Tune A is historically 
significant and makes use of Q^-ordered parton showers while the more recent tune PO makes 
use of pT-ordered showers. Both use the colour annealing model for colour reconnections. GAL 
is based on the tune SO [65] but uses the Generalised Area Law colour reconnection model. 
The different tunes show large differences for the considered observables which provide therefore 
important additional information to be taken into account in tunes. 

5.2. Double parton scattering in j + 3 jet events 

The double parton scattering fraction is measured [UB] in dependence of the second leading jet 
transverse momentum based on an integrated luminosity of 1.0 fb^^. The cross section for DPS 
is given by a^p = m ■ ^f^, where aA and as are the cross sections of the processes A and 
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Table 1. Parameters of different spatial parton density models calculated from aes- 



B, m is a constant permutation factor which corresponds to one in case of indistinguishable 
processes A and B. It assumes the value m = 2 in the case of distinguishable processes A and 
B as discussed here, aes is a process independent scaling parameter which gives a measure for 
the size of the effective interaction region. The expression as/^aes can be interpreted as the 
probability of a second interaction B given that a first one A has taken place. 

DPS is measured in 7 ± 3 jet events which benefits from the better energy measurement of 
the photon and a larger fraction of Double Parton (DP) events compared to 4-jet events [67j. 
Events with an isolated high transverse momentum photon of 60 < pT < 80 GeV are selected 
to increase the photon purity. Furthermore this requirement ensures a clean separation of the 
jet produced in the same parton scattering and jets produced by a different scattering. The 
jet transverse momenta are corrected to the hadronic final state. The first jet in transverse 
momentum is required to fulfil pT > 25 GeV and the second and third jet pT > 15 GeV. Jets 
are reconstructed with the iterative midpoint cone jet algorithm [3] with cone size R = 0.7. The 
main background is Single Parton (SP) scattering. Also Double Interaction (DI) events with two 
scatterings of two different pp collisions in the same bunch crossing (referred to as pile-up) have 
to be considered. They can be distinguished by the separate collision vertex. The contribution 
from single and double diffraction events represents < 1% of the total dijet cross section. The 
Pt spectrum for jets from dijet production falls faster than that for jets resulting from initial or 
final state radiation in the 7± jets events. Therefore DP fractions depend on the jet pT and the 
DP fractions as well as a^s can be determined in data by a set of equations, relating the different 
contributions in the three second jet transverse momentum bins 15 - 20, 20 - 25 and 25 - 30 GeV. 

In each momentum bin the discriminating variables Sp^, = "^y''"^^ + 

9 , - 1 J( \PTh,i)\ y , ( IPtUM Y . o _ 1 l( A<pi^,i) Y I f A^(j,fc) \2 p^hiated 
^p't - V2^[\p}\+\p^\) ^[\P^\+\P^\) ^"""^ "^-^ - v/2VU^(7,i)J +{s<p{j,k)) are evaluated, 

where (p is the azimuthal angle and 5 indicates the uncertainty of a given variable, defined 
between pairs of photon (7) and jet {i,j,k) objects. These variables are computed for the 
reconstructed particle pairing which minimises S. For this pairing and each variable the final 
observable AS = A(l){PT{'y,i), PtU, k)) is being determined. 

The measurement reveals that the DP fraction drops from 0.47 in the second jet 15 < 
Pt < 20 GeV transverse momentum bin to 0.23 in the second jet 25 < pT < 30 GeV 
transverse momentum bin. Averaging over the three second jet transverse momentum bins yields 
< CTeff >= 16.4 ± 0.3(stat) ± 2.3(syst) mb. Good agreement with previous CDF measurements 
in 4-jet [68| and 7±3 jet |69j [TO] events is observed. From the parameter a^s the proton radius 
in different parton spatial density models (see Table [1]) can be calculated. 
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5.3. Azimuthal decorrelations and double parton scattering in 7 + 2 jet and 'y + 3 jet events 
The azimuthal decorrelations and double parton scattering in 7 ± 2 jet and 7 ± 3 jet events are 
measured ^1] based on an integrated luminosity of 1.0 fb~^. This analysis provides an extension 
of the DPS analysis discussed above in allowing also for 7 ± 2 jet events, where a third jet is not 
reconstructed or did not pass the selection criteria. The gained statistics allows one to subdivide 
the cross section measurement in bins of the second jet transverse momentum and increases thus 
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Figure 6. Differential cross section in 7 + 2 jet events as a function of A(j) (left) and differential 
cross section in 7 + 3 jet events as a function of AS" (right) in the second leading jet bin 
respectively. The data is compared to various models. 



the sensitivity to MPI models. 

Events with an isolated high transverse momentum photon of 50 < pr < 90 GeV are selected. 
The jet transverse momenta are corrected to the hadronic final state. The first jet in transverse 
momentum is required to fulfil pT > 30 GeV and one or two more jets with > 15 GeV 
have to be found. Jets are reconstructed with the iterative midpoint cone jet algorithm ^ with 
cone size R = 0.7. The event selection for photons and jets is identical to the requirements 
of the previous analysis discussed above. In addition a separation requirement between any 
pair of selected objects among the photon and the jets of AR > 0.9 is used. The sample 
of DP candidates is selected from events with a single reconstructed pp collision vertex. For 
7 + 2 jet events the differential cross section as a function of the azimuthal angle difference 
between the photon + leading px jet system and the second leading pT jet A(/)(7 + jeti,jet2) is 
considered in three bins of second leading pT jet, namely 15 -20 (see Fig. [6l left), 20 - 25 and 
25 - 30 GeV. For 7 + 3 jet events the differential cross section as a function of the azimuthal 
angle difference between the photon + leading pj- jet system and the second + third leading 
Pt jet AS = A(/)(Pr(7, jet^), PT(jet2, jet3)), is measured only in a single 15-30 GeV second 
transverse momentum jet bin (see Fig. El right) for statistical reasons. The differential cross 
sections are normalised for an improved sensitivity to MPI models. Various models and tunes 
of PYTHIA ^60j and SHERPA [29j are compared to data. There are large differences between 
the models and the data confirm the presence of DP events. The data are close to Perugia 
(PO) [63], SO [65] and SHERPA [29] with MPI tunes, while predictions from previous PYTHIA 
MPI models with tunes A and DW, making use of Q^-ordered parton showers, are disfavoured. 

6. Conclusions 

A multitude of D0 measurements in the regime of perturbative and non-perturbative QCD has 
been elaborated and confronted to theory predictions. While pQCD predictions are in general in 
agreement with the various measurements, limitations of the non-perturbative parts (e.g. PDF's, 
photon fragmentation, hadronisation, underlying event) cannot be hidden and the experimental 
results can be exploited to optimise phenomenological models and their uncertainties. 
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